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Abstract 

In this note, we consider a nonlinear diffusion equation with a 
bistable reaction term arising in population dynamics. Given a rather 
general initial data, we investigate its behavior for small times as the 
reaction coefficient tends to infinity: we prove a generation of interface 
property. 
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1 Introduction 

We consider the degenerate parabolic problem 



' u t = A(u rn ) + \f(u) in Q T := £1 x (0, T) 
- = 



d(u r 
du 

u{x,Q) = u (x) 



on dn x (0, T) 
in Q , 



with e > a small parameter. Here Q is a smooth bounded domain in 
(N > 2), v is the Euclidian unit normal vector exterior to dVt and m > 2. 

We assume that / is smooth and has exactly three zeros < a < 1 such 
that 

/'(0)<0, /'(o)>0, /'(1)<0. (1.1) 
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Moreover we suppose that the initial function uq € C 2 (Q) is nonnegative, 
with support SuppuQ := {x £ fi| uo(x) > 0} CC 0. Furthermore we define 
the initial interface Tq by 

To := {x G fi| tio(x) = a} , 

and suppose that To is a smooth hypersurface without boundary such that 

T CC n and Vu (x) ^ if x G T , (1.2) 

uq > a in Oq 1 ^ , no < a in Oq ^ , (1-3) 

where Q,^ denotes the region enclosed by To and O °^ the region enclosed 
between d£l and IV 

We prove a generation of interface property, namely that the solution u £ 
quickly becomes close to 1 or 0, except in an 0{e) neighborhood of the initial 
interface To, creating a steep transition layer around IV More precisely, we 
are in presence of the following phenomenon: in the very early stage, the 
nonlinear diffusion term is negligible when compared with the reaction term 
e~ 2 f(u). Hence, under the rescaling in time r = t/e 2 , the equation is well 
approximated by the ordinary differential equation u T = f{u). In view of 
the bistable nature of /, u £ quickly approaches the stable equilibria of the 
ordinary differential equation, or 1, and an interface is formed between 
the regions {u £ ~ 0} and {u £ ~ 1}. 

The organization of this note is as follows. In Section [2] we briefly explain 
how Problem (P £ ) arises in population dynamics. In Section [31 we recall 
known results about the well-posedness of Problem (P £ ) and a compari- 
son principle. In Section [H we prove the generation of interface property 
for Problem (P £ ). To that purpose we construct sub- and super-solutions 
by modifying the solution of the corresponding ordinary differential equa- 
tion ut = e~ 2 f(u). We also show the optimality of the generation time 
t £ := /'( a ) _l£2 | l ne l an d prove that the thickness of the interface is of order 
0(e) at the generation time t £ . Our method of proof follows the same lines 
as that of [2] and p] . It is slightly different from those of Xinfu Chen [5] and 
[B], who transforms the reaction function /. We postpone to future work 
the study of the interface motion after the generation time of the interface. 

Finally let us mention articles involving the singular limit of reaction- 
diffusion equations with nonlinear diffusion. Feireisl [8] studies the singular 
limit of a degenerate parabolic equation in the whole space M N . He studies 
the problem in the scaling 

ut = eA(u m ) + -f(u) , (1.4) 
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where, in the limit e — > 0, the limit free boundary moves according to motion 
by constant speed (that of a related traveling wave). In a similar scaling, 
Hilhorst, Kersner, Logak and Mimura [TT] investigate the singular limit of 
this equation in a bounded domain with a monostable reaction term. In both 
of these papers, they prove that the solution u e of the nonlinear diffusion 
equation converges to or 1 on both sides of an interface moving with 
constant normal velocity. In this scaling, proofs about the interface motion 
can be performed with using only one term in the asymptotic expansion 
whereas we would need to use two terms in the case of Problem (P e ) as 
well as a suitable linearization procedure; this is far from trivial here since 
Problem (P £ ) is degenerate parabolic. 

2 The biological context 

In this section, we discuss nonlinear diffusion in population dynamics. It is 
well-known that the control of a population can be achieved by introducing 
density dependent birth or death rates. In [9], Gurney and Nisbet show that 
the introduction of a nonlinearity into the dispersal behavior of a species 
- which behaves in an otherwise linear way — can, in an inhomogeneous 
environment, lead to a regulatory effect. More precisely, they consider the 
equation 



where u(x,t) denotes the population density, G = G(x) the growth func- 
tion only depending on the location and t) the local population current 
density. By using the well-known random motion model one obtains the lin- 
ear equation ut = Au + G{x)u. Another possibility is to choose the biased 
random motion model where movements are largely random but slightly 
modified by the distribution of the fellows; the corresponding equation is 
then written as ut = Ait + div(itgradit) + G(x)u. Nevertheless, Carl 0] has 
observed that arctic ground squirrels migrate from densely populated areas 
into sparsely populated ones, even when the latter is less favorable (burrow 
sites not available, intensive predation). For such species, migration to avoid 
crowding, rather than random motion, is the primary cause of dispersal. To 
describe such movements, Gurney and Nisbet use the directed motion model 
where individuals can only stay put or move down the population gradient; 
this model yields the degenerate parabolic equation 



In [9], the authors perform a qualitative analysis of the three different dis- 
persal models (random motion, biased random motion and directed motion). 
They conclude that the introduction of density dependent diffusion is effi- 
cient to study the dynamics of a population which regulates its size below 



ut = — div j + G(x)u 




(2.1) 
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the carrying capacity set by the supply of nutrients. 

Gurtin and Mac Camy [10] proposed the class of equations which we 
study here and which involves degenerate diffusion and nonlinear reaction, 
namely 

u t = A{u m ) + /(it) , m>2. (2.2) 

In absence of a reaction term, equation (|2.2|) reduces to the so-called porous 
medium equation 

u t = A(u m ) , (2.3) 

which describes, among others, the flow of an ideal gas in a homogeneous 
medium {m > 2), groundwater infiltration (Boussinesq's equation, m = 2), 
the spread of a thin viscous film under gravity (m = 4), and thermal prop- 
agation in plasma (m ~ 6). The porous medium equation has been exten- 
sively investigated in literature: we refer to the book of Vasquez |12] and 
the references therein. The main feature of these equations is that they 
degenerate at the points where u = 0. As a consequence, a loss of regular- 
ity of solutions occurs and disturbances propagate with finite speed. This 
phenomenon contrasts with the infinite speed of propagation of solutions of 
the heat equation ut = Au. 



3 Comparison principle and well-posedness 

Since the diffusion term degenerates at the points where u e = 0, u £ is not 
smooth. This leads us to define a notion of weak solution for Problem (P e ), 
in a similar way as it is done by Aronson, Crandall and Peletier [3] for a 
corresponding one-dimensional problem. 

Definition 3.1. A function u £ : [0, oo) — > L 1 (J7) is a weak solution of 
Problem (P e ) if, for all T > 0, 

(i) u e G C ([0, oo); L 1 ^)) n L°°(Q r ) ; 

(ii) u e satisfies the integral equality 



u £ (TMT) - / / (u^ t + (n £ ) m A^) = / zw(O) 

Q, J JQt Jfl 



±f(u*)tp, (3.1) 



for all <f € C 2 (Qt) such that ip > and — — = on d£l. 

ov 

A sub-solution (respectively a super-solution) of Problem (P e ) is a function 
satisfying (i) and (ii) with equality replaced by < (respectively >). 
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Theorem 3.2 (Existence and comparison principle). Let T > be arbitrary. 
The following properties hold: 

(i) Let u~ and u + be a sub-solution and a super- solution of Problem (P £ ) 
with initial data and Uq respectively. 

If Uq < Uq then u~ < u + in Qt ; 

(ii) Problem (P £ ) has a unique weak solution u £ which is such that 

< u £ < max(l, HhoIIl 00 ^)) * n Qt (3-2) 

(iii) u £ € C(Qt~). 

The proof of Theorem 13.21 is standard; it can be performed by using the 
same lines as that of Theorem 5 in [3] . The continuity of u £ follows from [7] . 

The following result turns out to be an essential tool when constructing 
smooth sub- and super-solutions of Problem (P £ ). 

Lemma 3.3. Let u £ be a continuous nonnegative function in Q x [0, T]. 
Define n s t upp = {x € Q\ u £ {x,t) > 0} and T s t upp = dn s t upp for all t G [0,T]. 
Suppose that the family T := Uo < t<T^t UPP x W * s sufficiently smooth and 
let i/^ upp be the outward normal vector on Y^. upp . Suppose moreover that 

(i) V(u £ ) m is continuous in 0, x [0, T] ; 

(ii) C[u £ ] := uf - A(u £ ) m - \f{u £ ) = in {(x, t)eTlx [0, T] 

such that u £ (x, t) > 0} ; 

(iii) S£ = on dtf t upp , for all t e [0,T] . 

Then u £ is a solution of Problem (P £ ). Similarly u is a sub-solution (respec- 
tively a super- solution) of Problem (P £ ) if the equality in (ii) is replaced by 
< (resp. >) and if the equality in (iii) is replaced by < (resp. >). 

We refer to [11] for the proof. 

4 Generation of interface 

In this section we prove that, given a nearly arbitrary initial function uo, the 
solution u £ quickly becomes close to 1 or 0, except in an 0{e) neighborhood 
of the initial interface Tq, creating a steep transition layer around Tq. The 
time needed to develop such a transition layer is of order C(e 2 | lne|). 
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Theorem 4.1 (Generation of interface). Assume m > 2. Lei 7 G (0,min(a, 
1 — a)) be arbitrary and define fx as the derivative of f(u) at the unstable 
equilibrium u = a, namely 

H = f{a). (4.1) 

Moreover, set 

t e := n~ 1 e 2 \ lne| . 

Then there exist positive constants £0 and Mq such that, for all e G (0,£o), 

(i) for all x G Q, we have that 

< u £ (x,t £ ) < 1 + 7; 

(ii) for all x £ 0, such that \uq(x) — a\ > Mq£, we have that 

if u (x)>a + M e then u e (x, t e ) > 1 - 7 , (4.2) 
if u (x)<a-M e then u e (x,t e )<j. (4.3) 

Theorem 14.11 will be proved by constructing a suitable pair of sub- and 
super-solutions. As mentioned above, the nonlinear diffusion term is negli- 
gible in this early stage so that the behavior of the solution is governed by 

the ordinary differential equation ut = —^f(u). 

e z 

An immediate consequence of Theorem 14.11 is that u e (x,t £ ) is close to 
or 1, except in an 0(e) neighborhood of the initial interface Tq. In other 
words, the transition layers which have developed have an O(e) thickness. 

Corollary 4.2 (Thickness of the transition layers at time t e ). Let 77 G 
(0, min(o, 1— a)) be an arbitrary constant. Then there exist positive constants 
£0 and C such that, for all e S (0,£o), 

'[0,l + ? ? ] if xeMce(T ) 

u £ {x,t £ ) G \ [0,r/] if xettf\Nce{T Q ) (4.4) 

t [1-77, 1 + 7?] if xen^\M C e(r ), 

where denotes the region enclosed by Fq, Q,^ the region enclosed be- 
tween d£l and Tq, and 

M r (T ) := {x G f2, dist(x,T ) < r} 

denotes the r -neighborhood o/To- 

We will also show that the generation time t e := / u~ 1 £ 2 | lne\ is optimal. 
In other words, the interface is not fully developed until t becomes close to 
t £ . More precisely, the following result holds. 

Proposition 4.3. Denote by t^ in the smallest time such that (14. 4h holds. 
Then there exists a constant b = b(C) such that 

*min > M _1 e 2 (|lne| -b) 

for all £ G (0, £q). 
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4.1 Proof of the generation of interface property 
4.1.1 The bistable ordinary differential equation 

Let us first consider the problem without diffusion, namely 

ut = -of(u), u(x,0) = u (x) . 
Its solution can be written in the form 

u(x,t) = Y f-^, u (x) 

where V(t, £) is the solution of the ordinary differential equation 
f Y t (t,0 =f(Y(r,0) forr>0 



(4.5) 



Here £ ranges over the interval (—Co, Co), where Co := ||ito||z°°(n) + 1- We 
claim that Y has the following properties. 

Lemma 4.4. There exists a positive constant C such that the following holds 

(i) // f > then Y(t, £) > , 
// f < then Y(t, f ) < ; 

(ii) \Y(t,0\ <C ; 

(iii) y $ (r,0 >0; 

(iv) \^( T ,0\<C(e^-l), 

for all t > and aZZ £ S (—Co, Co). 

Properties (i) and (ii) are direct consequences of the profile of / - 
more precisely of the sign conditions / > in (— oo,0) U (a, 1) and / < 
in (0, a) U (1, oo) — and of the qualitative properties of the solution of the 
bistable ordinary differential equation (14.50 : for proofs of (iii) and (iv) we 
refer to [2], subsection 3.1. □ 

4.1.2 Construction of sub- and super-solutions 

We use the notation a + = max(a, 0). The sub- and super-solutions are given 
by 



wf(x,t) 



y(±, u (*)±£ 2 C*(e^/ £2 -l)) 



(4.6) 
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Lemma 4.5. There exist positive constants Eq and C* such that, for all e S 
(0, £q), (w~ , wf) is a pair of sub- and super- solutions for Problem (P £ ), in 
the domain $7 x [0,/x _1 e 2 | lne|]. Moreover, since also w~(x,0) = wf(x,0) = 
uq(x), it follows that 

w~(x,t) <u £ (x,t) <w+{x,t) for all (x,t) € tt x [0,^ _1 e 2 1 In e|] . (4.7) 

Proof. In order to prove that (w~ ,wf) is a pair of sub- and super-solutions 
for Problem (P £ ) for a suitable choice of £o and C g , we check that the 
sufficient conditions in Lemma 13.31 are satisfied. 

As for the sub-solution w~ , we remak that property (i) in Lemma 14.41 
implies that, for all t > 0, 

n s t upp [w-} = {xefl\ u (x) > e 2 C g {e^/ £2 - 1)} 

Tl upp [w-] := dn s t upp [w-} = {xen\ u (x) = e 2 C g {e^ £2 - 1)} . 

Choose (xo,io) such that xq € T^ pp [w~]; for (x,t) such that x € VL s t upp [w~\ 
we have 

V(w-) m (x,t) = mY m ~% (J^, u (x) - e 2 C*(e^ £2 - 1)) X7u (x) . 

Since Y(t, 0) = the equality above implies 

Urn V(w-r(x,t) = 0. 

(x,t)— >(xo,to) 
x€Q s t upp [we] 

Therefore conditions (i) and (iii) of Lemma 13.31 are satisfied by the sub- 
solution. 

As for the super-solution wf, we remark that property (i) of Lemma [4.41 
implies that, for all t > 0, 

nt up >t] = n 

rf pp [w+] -.= dn s t upp [ w +] = on . 

Hence condition (iii) of Lemma 13.31 for the super-solution is a direct conse- 
quence of the fact that SuppuQ CC O, whereas condition (i) is obviously 
satisfied. 

It remains to prove that 

C[w-] := (w-) t - A(w-) m ~ ^f{w~) < 0, 

in {(x,t) £ Q, x [0, fi~ 1 e 2 \lne\] such that w~(x,t) > 0} and that C[wf] > 
in {(x,t) eTlx [O^-^^lnel]}. 
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In view of the ordinary differential equation (14, 5D . straightforward com- 
putations yield 



-Yp 



C* /ie^ 2 + m{m - l)Y m - 2 Y^\Vu \ 

Ye. 



+ mY m ~ 1 ^\Vu \ 2 + mY m - 1 Au , (4.8) 



Yt 



in Q s t upp [w £ ], where the function Y and its derivatives are taken at the 
point (r, £) = (t/e 2 ,uo(x) — e 2 C*(e M< / e — 1)). Moreover since the term 

n 



m(m — l)Y m 2 l^ 2 |Vuo| 2 is nonnegative, it follows that 



C* ue^ t/e2 +mY m ~ 1 ^- \Vu \ 2 + mY^A 



Yt 



UQ 



We note that, in the range < t < /Lt _1 e 2 |lne|, we have, for eq sufficiently 
small, 

£ = u (x) - e 2 C*(e^ 2 - 1) e (-Co, C ) • 

We deduce from the properties (ii)-(iv) stated in Lemma [4.41 that there exist 
positive constants C\ and C*2 — only depending on m, Co, C, ||V«o||L°°(f2) 
and || Aito 1 1 loo — such that 



£[w-} < -y e 



(C*n - Ci)e^/ £2 - C 2 



which implies that C[w £ ] < if C* is chosen large enough. 
As for the super-solution we obtain 



Y f 



C* ne^ t/£2 - m(m - l)y m " 2 y, |Vu | 2 



m ym-l IT7„._ |2 ™ \rm-l , 



Yt 



[V«o - mY m ~ L Au Q 



, (4.9) 



and the assumption that m > 2 gives an upper bound for |Y m_2 |. Following 
the same argument as above one can prove that > for C* sufficiently 

large. This completes the proof of Lemma 14.51 □ 

4.1.3 Proof of Theorem 14.11 

In order to prove Theorem l4.1l we first present basic estimates of the function 
Y after a time of order r ~ | lne|. 

Lemma 4.6. Let 7 £ (0, min(a, 1 — a)) be arbitrary. There exist positive 
constants Eq and Cy such that, for all e € (0,£o), 

(i) for all £e (-Co, Co), 



7 <y(^" 1 |lneU) < 1 + 7; 



(4.10) 
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(ii) for all £ € (—Co, Co) such that |£ — a\ > Cye, we have that 

if ^>a + C Y e then Y{^~ 1 \ lne|,£) > 1 - 7 (4.11) 
*/ t<a-C Y e then Y{^~ 1 j lne|, £) < 7 . (4.12) 

These estimates illustrate the stability of the equilibria and 1 for the 
bistable ordinary differential equation (|4.5p . For more details we refer the 
reader to the proof of Lemma 3.9 in [2]. □ 

We are now ready to prove Theorem 14.11 By setting t = fi~ 1 e 2 \ lne| in 
(|4.7p . we obtain 

[Y (^ _1 |hie|,uo(ac) - (C*e - C*e 2 ))] + 

< u £ {x,n~ l e 2 \lne\) < [Y (/x" 1 ] lne\,u {x) + C*e - C*e 2 )] + . (4.13) 

Since, for eo small enough, uq(x) + (C*e — C*e 2 ) € (0, Co), the assertion (i) 
of Theorem 14.11 is a direct consequence of (14.101) and (|4.13j) . 

Next we prove (|4,2p . We choose Mq large enough so that Mqe — C*e + 
C*e 2 > Cye. Then, for all x € Q such that uo(x) > a + Mo£, we have that 
ii (a;) - (C*e - C*e 2 ) > a + Cye, which we combine with (|4~T3j) and ([17TT]) 
to deduce that 

^(a;,//" 1 ^ 2 ! lnej) > 1 - 7. 

The inequality (14. 3p can be shown in a similar way. This completes the proof 
of Theorem 14.11 □ 

Remark 4.7. Theorem 14.11 remains true if we perturb the reaction function 
f(u) by order e, setting for instance f{u) = f(u) — eg(x,t,u). To deal 
with this more general case, we proceed as follows. We first consider a 
slightly perturbed reaction function, namely fs(u) = f(u) + 5 which, for 5 
small enough, is still of bistable type. Define a(5) as its unstable zero and 
/j,(6) := f'(a(5)) as the slope of f$ in this point. We then define Y(t, £; 6) as 
the solution of the initial value problem 

f Y T (r, £ 6) = f s (Y(r, f ; 5)) for r > 



Finally we construct a pair of sub- and super-solutions in the form 

w±(x,t) 



Y I -4, uo(x) ± £ 2 r(±e£, ^); ±e£ 



where r(<5, r) := C*(e^^ T — 1) and C? := H^Hi 00 (o,Cb) - ^or more details 
and proofs we refer to [2], more precisely to Section 4 which deals with a 
generation of interface property for an equation with linear diffusion and an 
unbalanced reaction term f(u). □ 
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4.2 Proof of the optimality of the generation time 

We show below that the generation time t £ := fj, e 2 \hie\ is optimal. In 
other words, the interface is not fully developed before t is close to t £ . We 
will need the following lemma about the solution Y(t, £) of the corresponding 
bistable ordinary differential equation. 

Lemma 4.8. Let ij € (0, min(a, 1 — a)) be arbitrary. Then there exist positive 
constants C\ = C\{rf) and C2 = C^ 7 ?) such that, if £ € (a, 1 — rf) then, for 
every r > such that Y{r,£f) remains in the interval (a, 1 — rf), we have 

de^iC - a) < Y(t,0 - a < C 2 e' tT '(£ - a) . (4.15) 

We refer to [2] Corollary 3.5 for the proof of Lemma 14.81 □ 

Proof of Proposition 14.31 For each b > 0, we set 

t £ (b) := fi- 1 e 2 (\lns\ -6), 

and evaluate u £ (x,t £ (b)) at a point x E SIq such that <iisi(2;, To) = Ce. 
Define C := H^oHc 2 ^)- Since uq = a on To, we have that 

«o(a?) < a + CCe, (4.16) 
which implies together with Lemma 14.81 that 

wf (x, t e (6) ) = y (V 1 ( I In e \ - b) , u (x) + eC*e~ b - e 2 C*) 

< a + C 2 el ln£ l- fe (u (x) + eC*e" fe - e 2 C* - a) 

< a + C 2 e- 1 e- 6 (CCe + eC* e - b ) 
= a + Cae-^CC + ^e" 6 ). 

Now choose b large enough so that 

a + C 2 e- b (CC + C*e' b ) < 1 -77; (4.17) 

the inequalities (|4.17p and (|4.7p then yield 

u £ (x,t £ (b)) < wf(x,t £ (b)) <l-7]. 

Therefore (|4.4|) does not hold at i = t £ (b), and hence t £ (fe) < t^ in . This 
completes the proof of Proposition 14.31 □ 
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